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Abstract. We study the C ARDPath constraint. This ensures a given 
constraint holds a number of times down a sequence of variables. 
We show that SLIDE, a special case of CardPath where the slid 
constraint must hold always, can be used to encode a wide range 
of sliding sequence constraints including CardPath itself. We con- 
sider how to propagate SLIDE and provide a complete propagator for 
CardPath. Since propagation is NP-hard in general, we identify 
special cases where propagation takes polynomial time. Our experi- 
ments demonstrate that using SLIDE to encode global constraints can 
be as efficient and effective as specialised propagators. 



1 INTRODUCTION 



In many scheduling problems, we have a sequence of decision vari- 
ables and a constraint which applies down the sequence. For exam- 
ple, in the car sequencing problem, we need to decide the sequence 
of cars on a production line. We might have a constraint on how often 
a particular option is met (e.g. 1 out of 3 cars can have a sun-roof). As 
a second example, in a nurse rostering problem, we need to decide 
the sequence of shifts worked by nurses. We might have a constraint 
on how many consecutive night shifts any nurse can work. Such con- 
straints have been classified as sliding sequence constraints (7). To 
model such constraints, we can use the CardPath constraint. This 
ensures that a given constraint holds a number of times down a se- 
quence of variables |5|. We identify a special case of CardPath 
which we call SLIDE, that is interesting for several reasons. First, 
many sliding sequence constraints, including CardPath, can easily 
be encoded using this special case. SLIDE is therefore a "general- 
purpose" constraint for encoding sliding sequencing constraints. This 
is an especially easy way to provide propagators for such global con- 
straints within a constraint toolkit. Second, we give a propagator for 
enforcing generalised arc-consistency on SLIDE. By comparison, the 
previous propagator for CardPath given in |5| does not prune all 
possible values. Third, SLIDE can be as efficient and effective as spe- 
cialised propagators in solving sequencing problems. 
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2 CARDPATH AND SLIDE CONSTRAINTS 

A constraint satisfaction problem consists of a set of variables, each 
with a finite domain of values, and a set of constraints specifying 
allowed combinations of values for given sets of variables. We use 
capital letters for variables (e.g. X), and lower case for values (e.g. 
d). We write D{X) for the domain of variable X. Constraint solvers 
typically explore partial assignments enforcing a local consistency 
property. A constraint is generalised arc consistent (GAC) iff when 
a variable is assigned any value in its domain, there exist compatible 
values in the domains of all the other variables of the constraint. 

The CardPath constraint was introduced in |5|. If C is a con- 
straint of arity k then CardPath(A'', [Xi,. . . C) holds iff 
C{Xi, . . . , Xi+t-i) holds A'' times for l<i<n — fc + 1. For 
example, we can count the number of changes in the type of shift 
with CARDPATH(Af, [Xi, . . . , Xn], Note that CardPath can 
be used to encode a range of Boolean connectives since A'^ > 1 
gives disjunction, A'^ = 1 gives exclusive or, and N = gives nega- 
tion. We shall focus on a special case of the CardPath constraint 
where the slid constraint holds always. Slide(C, [Xi, . . . , X„]) 
holds iff C{Xi, . . . ,Xi+fc_i) holds for all 1 < i < n - fc + 1. 
That is, a CardPath constraint in which A'^ = n — + We 
also consider a more complex form of SLIDE that applies only ev- 
ery j variables. More precisely, SLIDEj(C, [Xi, . . . , Xn]) holds iff 
C{Xij+i,...,Xij+k) holds for < i < By definition 
SlidEj for j = 1 is equivalent to SLIDE. 

Beldiceanu and Carlsson have shown that CardPath can encode 
a wide range of constraints like CHANGE, SMOOTH, AmongSeq 
and SlidingSum |5 |. As we discuss later, SLIDE provides a simple 
way to encode such sliding sequencing constraints. It can also en- 
code many other more complex sliding sequencing constraints like 
Regular |T6|, Stretch | 13 |, and Lex [TJ, as well as many types 
of chanelling constraints like ELEMENT II9I and optimisation con- 
straints like the soft forms of REGtJLAR |20l . More interestingly, 
CardPath can itself be encoded into a SLIDE constraint. In [5), a 
propagator for CardPath is proposed that greedily constructs up- 
per and lower bounds on the number of (un)satisfied constraints by 
posting and retracting (the negation of) each of the constraints. This 
propagator does not achieve GAC. We propose here a complete prop- 
agator for enforcing GAC on SLIDE. SLIDE thus provides a GAC 
propagator for C ARDPath. In addition, SLIDE provides a GAC prop- 
agator for any of the other global constraints it can encode. As our 
experimental results reveal, SLIDE can be as efficient and effective 
as specialised propagators. 

We illustrate the usefulness of SLIDE with the AmongSeq 
constraint which ensures that values occur with some given fre- 
quency. For instance, we might want that no more than 3 out 
of every sequence of 7 shift variables are a "night shift". More 



precisely, AmongSeq(Z, u, k, [Xi, . . . , X„],v) holds iff between 
I and It variables in every sequence of k variables take value 
in the ground set v jS). We can encode this using SLIDE. More 
precisely, AMONGSBQ{l,u,k,[Xi, . . . ,Xn],v) can be encoded 
as SLIDE(Z)f [Xi,.. . ,X„]) where Df'^ is an instance of the 
Among constraint i;^. Di'^{Xi, . . . , Xi+fc_i) holds iff between / 
and u variables take values in the set v. For example, suppose 2 of 
every 3 variables along a sequence Xi . . . X5 should take the value 
a, where Xi = a and X2, ■ ■ ■ , X^ £ {a, fe}. This can be encoded 
as SLIDE(i;, [Xi, X2, X3, X4, X5I) where E{X,,X,+x,X,+-z) en- 
sures two of its three variables take a. This SLIDE constraint ensures 
thatS(Xi,X2,X3), E(X2,Xi,X^) andS(X3,X4,X5) all hold. 
Note that each ternary constraint is GAC. However, enforcing GAC 
on the Slide constraint sets X4 = a as there are only two satisfying 
assignments and neither have Xi^ = b. 



3 SLIDE WITH MULTIPLE SEQUENCES 

We often wish to slide a constraint down two or more sequences of 
variables at once. For example, suppose we want to ensure that two 
vectors of variables, Xi to Xn and Yi to Y„ differ at every index. We 
can encode such a constraint by interleaving the two sequences and 
sliding a constraint down the single sequence with a suitable offset. 
In our example, we simply post SlidE2(7^, [Xi, Yi, . . . , X„, Yn]). 
As a second example of sliding down multiple sequences of vari- 
ables, consider the constraint Regular(^, [Xi, . . . , Xn]). This en- 
sures that the values taken by a sequence of variables form a string 
accepted by a deterministic finite automaton ^ 1 1 6 1 . This global con- 
straint is useful in scheduling, rostering and sequencing problems to 
ensure certain patterns do (or do not) occur over time. It can be used 
to encode a wide range of other global constraints including: AMONG 
(D, Contiguity ||15J, Lex and Precedence fT4l . 

To encode the REGULAR constraint with SLIDE, we intro- 
duce variables, Qi to record the state of the automaton. We 
then post SlidE2(F, [Qo, Xi, Qi, . . . , X„, Q„]) where Qo is set 
to the starting state, Qn is restricted to accepting states, and 
F{Qi,Xi+i,Qi+i) holds iff Qi+i = S{Xi,Qi) where S is the 
transition function of the automaton. If we decompose this encoding 
into the conjunction of slid constraints, we get a set of constraints 
similar to |6|. Enforcing GAC on this encoding ensures GAC on 
Regular and, by exploiting functionally of F, takes 0(ndq) time 
where d is the number of values for Xi and q is the number of states 
of the automaton. This is asymptotically identical to the specialised 
Regular propagator 1 16|. This encoding is highly competitive in 
practice with the specialized propagator | 2]. 

One advantage of this encoding is that it gives explicit access to the 
states of the automaton. Consider, for example, a rostering problem 
where workers are allowed to work for up to three consecutive shifts. 
This can be specified with a simple REGULAR constraint. Suppose 
now we want to minimise the number of times a worker has to work 
for three consecutive shifts. To encode this, we can post an AMONG 
constraint on the state variables to count the number of times we 
visit the state representing three consecutive shifts, and minimise 
the value taken by this variable. As we shall see later in the exper- 
iments, the encoding also gives an efficient incremental propagator. 
In fact, the complexity of repeatedly enforcing GAC on this encoding 
of the Regular constraint down the whole branch of a backtracking 
search tree is just 0{ndq) time. 



4 SLIDE WITH COUNTERS 

We may want to slide a constraint on a sequence of variables comput- 
ing a count. We can use SLIDE to encode such constraints by incre- 
mentally computing the count in an additional sequence of variables. 
Consider, for example, CardPath(A'', [Xi, . . . , X„], C). For sim- 
plicity, we consider k — 2 (i.e., C is binary). The generalisation to 
other k is straightforward. We introduce a sequence of integer vari- 
ables Mi in which to accumulate the count. We encode CardPath 
as SlidE2(G, [Ml, Xi, . . . , M„,X„]) where Mi = 0, M„ = N, 
and G(Mi,Xi,Mi+i,Xi+i) is defined as: if C(Xi,Xi+i) holds 
then Mi+i = Mi + 1, otherwise Mi+i = Mi. GAC on SLIDE 
ensures GAC on CardPath. 

As a second example, consider the STRETCH constraint [131 . 
Given variables Xi to X„ taking values from a set of shift 
types r, a set n of ordered pairs from r x r, and functions 
shortest{t) and longest{t) giving the minimum and maximum 
length of a stretch of type t, Stretch([Xi, . . . , X„]) holds iff each 
stretch of type t has length between shortest{t) and longest{t); 
and consecutive types of stretches are in tt. We can encode 
STRETCHas SlidE2(H, [Xi, Qi, . . . , X„, Q„]) where Qi = 1 and 
f/(X„ X,+i,0„Q,+i) holds iff (1) Xi =X,+i,Q,+i = 1 + Q„ 
and Qi+i < longest{Xi); or (2) Xi / Xi+i, (Xi,Xi+i) £ tt, 
Qi > shortest(Xi) and Qi+i = 1. GAC on SLIDE ensures GAC 
on Stretch. 

5 OTHER EXAMPLES OF SLIDE 

There are many other examples of global constraints which we can 
encode using SLIDE. For example, we can encode LEX (7) us- 
ing Slide. Lex holds iff a vector of variables [Xi..X„] is lexico- 
graphically smaller than another vector of variables [Yi..F„]. We 
introduce a sequence of Boolean variables Bi to indicate if the 
vectors have been ordered by position i — 1. Hence Bi — 0. 
We then encode LEX as SlidE3(J, [Bi,Xi, Yi, . . . ,Bn,Xn,Y„]) 
where I{B,,X„ Y,B,+i) holds iff (B, = B,+i = A X, = Yi) 
or (B, = A B^+i = 1 A X, < Y,) or [B, = B,+ i = 1). 
This gives us a linear time propagator as efficient and incremen- 
tal as the specialised algorithm in |I2|. As a second example, we 
can encode many types of channelling constraints using SLIDE like 
Domain (TT], LinkSet2Booleans (T) and Element (i3- As 
a final example, we can encode "optimisation" constraints like the 
soft form of the REGULAR constraint which measures the Hamming 
or edit distance to a regular string |20|. There are, however, con- 
straints that can be encoded using SLIDE which do not give as ef- 
ficient and effective propagators as specialised algorithms (e.g. the 
global AllDifferent constraint |18 |). 

6 PROPAGATING SLIDE 

A constraint like SLIDE is only really useful if we can propagate 
it efficiently and effectively. The simplest possible way to propa- 
gate SLIDEj(C, [Xi, . . . , X„]) is to decompose it into a sequence 
of constraints, C(Xij+i, . . . , X,j+k) for < i < ^ and let the 
constraint solver propagate the decomposition. Surprisingly, this is 
enough to achieve GAC in many cases. For example, we can achieve 
GAC in this way on the SLIDE encoding of the REGULAR constraint. 
If the constraints in the decomposition overlap on just one variable 
then the constraint graph is Berge acyclic |4|, and enforcing GAC on 
the decomposition of SLIDEj achieves GAC on SlidEj. Similarly, 
enforcing GAC on the decomposition achieves GAC on SLIDEj if 



the constraint being slid is monotone. A constraint C is monotone 
iff tliere exists a total ordering -< of the values such that for any two 
values v,w, if V -< w then v can replace w in any support for C. For 
instance, the constraints AMONG and SUM are monotone if either no 
upper bound, or no lower bound is given. 

Theorem 1 Enforcing GAC over each constraint in the decomposi- 
tion o/SlidEj achieves GAC on SLIDEj if the constraint being slid 
is monotone. 

Proof: For an arbitrary value v £ D(X), we show that if every 
constraint is GAC, then we can build a support for X = v on SLIDEj . 
For any variable other than X, we choose the smallest value in the 
total order. This is the value that can be substituted for any other value 
in the same domain. A tuple built this way satisfies all the constraints 
being slid since we know that there exists a support for each (they are 
GAC), and the values we chose can be substituted for this support. □ 
In the general case, when constraints overlap on more than one 
variable (e.g. in the SLIDE encoding of AmongSeq), we need to do 
more work to achieve GAC. We distinguish two cases: when the arity 
of the constraint being slid is not fixed, and when the arity is fixed. 
We show that enforcing GAC in the former case is NP-hard. 

Theorem 2 Enforcing GAC on Slide(C, [Xi , . . . , X„]) is NP- 
hard when the arity of C is not fixed even if enforcing GAC on C 
is itself polynomial. 

Proof: We give a reduction from 3-SAT in A'' variables and M 
clauses. We introduce variables Xl for 1 < i < A'' + 1 and 
1 < J < For each clause j, if the clause is Xa V -^Xb V Xc, then 
we set XI £ {xa, ~'Xh, Xc\ to represent the values that make this 
clause true. For each clause j, we set Xf^-^ £ {0, 1} for 1 < i < A'^ 
to represent a truth assignment. Hence, we duplicate the truth 
assignment for each clause. We now build the following constraint 
Slide(C, [X-i X]^^-^, ..jXf, ..,Xlj^-^, ..,Xi , .., Xjv+i]) 
where C has arity A*' + 1. We construct C(Yi, . . . , Vjv+i) to hold 
iff Yi = Xd and Yi+d = 1, or Yi = -^Xd and Yi+d = 0. (in 
these two cases, the value assigned to Yi represents the literal that 
makes clause j true), or Yi G {0, 1} and Yi — Yi+M+i (in this 
case, the truth assignment is passed down the sequence). Enforcing 
GAC on C is polynomial and an assignment satisfying the Slide 
constraint corresponds to a satisfying assignment for the original 
3-SAT problem. □ 

When the arity of the constraint being slid is not great, we can en- 
force GAC on Slide using dynamic programming (DP) in a similar 
way to the DP-based propagators for the REGULAR and STRETCH 
constraints I16II13I . A much simpler method, however, which is just 
as efficient and effective as dynamic programming is to exploit a 
variation of the dual encoding into binary constraints IITOj based on 
tuples of support. Such an encoding was proposed in |T| for a par- 
ticular sliding constraint. Here we show that this method is more 
general and can be used for arbitrary SLIDE constraints. Using such 
an encoding, SLIDE can be easily added to any constraint solver. We 
illustrate the intersection encoding by means of an example. 

Consider again the AmongSeq example in which 2 of ev- 
ery 3 variables of Xi . . . X5 should take the value a, where 
Xl = a and X2,...,X5 £ {o.,b}. We can encode this as 
Slide(£, [Xi,X2,X3,X4,X5]) where E{X„ X,+ i, X,+2) is an 
instance of the AMONG constraint that ensures two of its three vari- 
ables take a. If the sliding constraint has arity k, we introduce an 
intersection variable for each subsequence of fc — 1 variables of 
Slide. The first intersection variable Vi has a domain containing 
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Figure 1. Intersection encoding 

all tuples from D{Xi) x . . . x D{Xk-i). The jth intersection vari- 
able Vj has domain containing D{Xj) x . . . x D(Xj+fc_2). And 
so on until Vn-k+2. In our example in Fig[T] this gives D{V\) = 
D{Xi) X D{X2), D{V4) = DiXi) X ^(Xs). We then post bi- 
nary compatibility constraints between consecutive intersection vari- 
ables. These constraints ensure that the two intersection variables 
assign {k — 1) -tuples that agree on the values of their k — 2 com- 
mon original variables (like constraints in the dual encoding). They 
also ensure that the fc-tuple formed by the two (fc — l)-tuples sat- 
isfies the corresponding instance of the slid constraint. For instance, 
in Fig[T] the binary constraint between Vi and V2 does not allow 
the pair {ab, aa) because the second argument of ah for Vi (value h 
for X2) is in conflict with the first argument of aa for V2 (value a for 
X2). That same constraint between Vi and V2 does not allow the pair 
{ab, bb) because the tuple abb is not allowed by E{Xi,X2, Xz). 

Enforcing AC on such compatibility constraints prunes aa and bb 
from V2, ab and bb from V3, and ba and bb from V4. Finally, we 
post binary channelling constraints to link the tuples to the original 
variables. One such constraint for each original variable is sufficient. 
For example, we can have a channelling constraint between V4 and 
Xi which ensures that the first argument of the tuple assigned to V4 
equals the value assigned to X4,. Enforcing AC on this channelling 
constraint prunes b from the domain of X4. We could instead post a 
channelling constraint between V3, and X4 ensuring that the second 
argument in V3 equals X4. The AmongSeq constraint is now GAC. 

Theorem 3 Enforcing AC on the intersection encoding of SLIDE 
achieves GAC in 0{nd!') time and 0{nd''~^) space where k is the 
arity of the constraint to slide and d is the maximum domain size. 

Proof: The constraint graph associated with the intersection 
encoding is a tree. Enforcing AC on this therefore achieves GAC. 
Enforcing AC on the channelling constraints then ensures that 
the domains of the original variables are pruned appropriately. As 
we introduce 0(n) intersection variables, and each can contain 
0{d^^^) tuples, the intersection encoding requires 0{nd''^^) 
space. Enforcing AC on a compatibility constraint between two 
intersection variables Vi and Vi+i takes 0{d'') time as each tuple 
in the intersection variable Vi has at most d supports which are the 
tuples of Vi+i that are equal to Vi on their fc — 2 common arguments. 
Enforcing AC on 0{n) such constraints therefore takes Oind^) 
time. Finally, enforcing AC on each of the 0(n) channelling 
constraints takes 0{df'~^) time as they are functional. Hence, the 
total time complexity is 0{nd''). □ 

Arc consistency on the intersection encoding simulates pairwise 
consistency on the decomposition. It does this efficiently as inter- 
section variables represent in extension 'only' the intersections. This 
is sufficient because the constraint graph is acyclic. This encoding 
is also very easy to implement in any constraint solver. It has good 
incremental properties. Only those constraints associated with a vari- 



able which changes need to wake up. 

The intersection encoding of SLIDEj for j > 1 is less expensive 
to build than for j = 1 as we need intersection variables for subse- 
quences of less than fc — 1 variables. For 1 < j < k/2, we introduce 
intersection variables for subsequences of variables of length k — j 
starting at indices l,j + 1, 2j + 1... whose domains contain {k — j)- 
tuples of assignments. Compatibility and channelling constraints are 
defined as with j = 1. If j > k/2, two consecutive intersection vari- 
ables (for two subsequences of fc — j variables) involve less than k 
variables of the SLIDEj. The compatibility constraint between them 
cannot thus ensure the satisfaction of the slid constraint. We therefore 
introduce intersection variables for subsequences of length \k/2~\ 
starting at indices l,j + l,2j + 1... and for subsequences of length 
[fe/2] finishing at indices fc, j + k, 2j + k... The compatibility con- 
straint between two consecutive intersection variables representing 
the subsequence starting at index pj + 1 and the subsequence fin- 
ishing at index pj + k ensures satisfaction of the {p+ l)th instance 
of the slid constraint. The compatibility constraint between two con- 
secutive intersection variables representing subsequence finishing at 
index pj + k and the subsequence starting at index (p + l)j + 1 
ensures the consistency of the arguments in the intersection of two 
instances of the slid constraint. 



choices do not affect results. We schedule the days in chronological 
order and within each day we allocate a shift to every nurse in lexi- 
cographical order. Initial experiments show that this is more efficient 
than the minimum domain heuristic. However, it restricts the variety 
of domains passed to the propagators, and thus hinders any demon- 
stration of differences in pruning. We therefore also use a more ran- 
dom heuristic. We allocate within each day a shift to every nurse 
randomly with 20% frequency and lexicographically otherwise. 
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7 EXPERIMENTS 

We now demonstrate the practical value of SLIDE. Due to space 
limits, we only report detailed results on a nurse scheduling prob- 
lem, and summarise the results on balanced incomplete block design 
generation and car sequencing problems. Experiments are performed 
with ILOG Solver 6.2 on a 2.8GHz Intel computer running Linux. 

We consider a Nurse Scheduling Problem (9) in which we gener- 
ate a schedule of shift duties for a short-term planning period. There 
are three types of shifts (day, evening, and night). We ensure that 

(1) each nurse takes a day off or is assigned to an available shift; 

(2) each shift has a minimum required number of nurses; (3) each 
nurse's work load is between specific lower and upper bounds; (4) 
each nurse works at most 5 consecutive days; (5) each nurse has at 
least 12 hours of break between two shifts; (6) the shift assigned to 
a nurse does not change more than once every three days. We con- 
struct four different models, all with variables indicating what type 
of shift, if any, each nurse is working on each day. We break symme- 
try between the nurses with lex concstraints. The constraints (l)-(3) 
are enforced using global cardinality constraints. Constraints (4), (5) 
and (6) form sequences of respectively 6-ary, binary and ternary con- 
straints. Since (4) is monotone, we simply post the decomposition in 
the first three models. This achieves GAC by Theorem[T] The mod- 
els differ in how (5) and (6) are propagated. In decomp, they are 
decomposed into conjunction of slid constraints. In amongseq, (5) 
is decomposed and (6) is enforced using the AmongSeq constraint 
of ILOG Solver (called IloSequence). The combination of (5) 
and (6) are enforced by SLIDE in slide. Finally, in slidec, we 
use Slide for the combination of (4), (5), and (6). 

We test the models using the instances available at 
^ttp.VAvww.projectmanageinent. iigent. be/nsp.php in which nurses 
have no maximum workload, but a set of preferences to optimise. 
We ignore these preferences and post a constraint bounding the 
maximum workload to at most 5 day shifts, 4 evening shifts and 2 
night shifts per nurse and per week. Similarly, each nurse must have 
at least 2 rest days per week. We solve three samples of instances 
involving 25, 30 and 60 nurses to schedule over 28 days. 

We use the same variable ordering for all models so that heuristic 



#solved bts time bts time 

25 nurses, 28 days (99 instances) 

decomp g5 — 35ngT 739 — mm nn35" 

amongseq 85 35401 14.43 35401 14.43 

slide 97 1699 1.00 1547 0.92 

slidec 97 457 0.58 438 0.56 

30 nurses, 28 days (99 instances) 
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60 nurses, 28 days (100 instances) 

decomp 1 11^55 Tim imn 142.90 

amongseq 2 122406 119.40 122406 119.40 

slide 27 562 0.65 2367 2.19 

slidec 34 542 3.96 1368 , 6.38 

Tabled. — Nur.>ie scheduling wllh random variable ordering (^ on Instances 

solved by all methods, ^ on instances solved by the method). 

Tables [T] and |2] report the mean runtime and fails to solve the in- 
stances with 5 minutes cutoff. Between the first three models, the best 
results are due to slide. We solve more instances with slide, 
as well as explore a smaller tree. By developing a propagator for 
a generic constraint like SLIDE, we can increase pruning without 
hurting efficiency. Note that slide always performs better than 
amongseq. A possible reason is that AmongSeq cannot encode 
constraint (6) as directly as SLIDE. As in previous models, we need 
to channel into Boolean variables and post AmongSeq on them. 
This may not give as effective and efficient pruning. SLIDE thus of- 
fers both modelling and solving advantages over existing sequencing 
constraints. Note also that slidec solves additional instances in the 
time limit. This is not suprising as the model slides the combination 
of the constraints (4), (5), and (6). Recall that the sliding constraint 
of (4) is 6-ary. It is pleasing to note that the intersection encoding 
performs well even in the presence of such a high arity constraint. 

We also ran experiments on Balanced Incomplete Block Designs 
(BIBDs) and car sequencing. For BIBD, we use the model in 1121 
which contains LEX constraints. We propagate these either using 
the specialised algorithm of 1 12| or the SLIDE encoding. As both 
propagators maintain GAC, we only compare runtimes. Results on 
large instances show that the SLIDE model is as efficient as the LEX 
model. For car sequencing, we test the scalability of SLIDE on large 
arity constraints and large domains using 80 instances from CSPLib. 
Unlike a model using IloSequence, our SLIDE model does not 
combine reasoning about overall cardinality of a configuration with 



the sequence of AMONG constraints. Hence, it is not as efficient: 26 
instances were solved with SLIDE within the five minute cutoff, com- 
pared to 39 with IloSequence. However, 9 of the instances solved 
with Slide were not solved by IloSequence. The memory over- 
head of the Slide propagator was not excessive despite the slid con- 
straints having arity 5 and domains of size 30. The SLIDE model used 
on average 22Mb of space, compared to 5Mb for IloSequence. 

8 RELATED WORK 

Pesant introduced the REGULAR constraint, and gave a propagator 
based on dynamic programming to enforce GAC [16]. As we saw, 
the Regular constraint can be encoded using a simple SLIDE con- 
straint. In this simple case, the dynamic programming machinery of 
Pesant's propagator is unnecessary as the decomposition into ternary 
constraints does not hinder propagation. We have found that SLIDE 
is as efficient as REGULAR in practice |2|. Furthermore, our encod- 
ing introduces variables for representing the states. Access to the 
state variables may be useful (e.g. for expressing objective func- 
tions). Although an objective function can be represented with the 
CostRegular constraint |TT], this is limited to the sum of the 
variable-value assignment costs. Our encoding is more flexible, al- 
lowing different objective functions like the min function used in the 
example in Section[3] 

Beldiceanu, Carlsson, Debruyne and Petit have proposed specify- 
ing global constraints by means of deterministic finite automata aug- 
mented with counters f6l . They automatically construct propagators 
for such automata by decomposing the specification into a sequence 
of signature and transition constraints. This gives an encoding sim- 
ilar to our Slide encoding of the REGULAR constraint. There are, 
however, a number of advantages of SLIDE over using an automaton. 
If the automaton uses counters, pairwise consistency is needed to 
guarantee GAC (and most constraint toolkits do not support pairwise 
consistency). We can encode such automata using a SLIDE where we 
introduce an additional sequence of variables for each counter. SLIDE 
thus provides a GAC propagator for such automata. Moreover, SLIDE 
has a better complexity than a brute-force pairwise consistency algo- 
rithm based on the dual encoding as it considers only the intersection 
variables, reducing the space complexity by a factor of d. 

Hellsten, Pesant and van Beek developed a GAC propagator for 
the Stretch constraint based on dynamic programming similar to 
that for the REGULAR constraint [ 13 1. As we have shown, we can en- 
code the Stretch constraint and maintain GAC using SLIDE. Sev- 
eral propagators for the AmongSeq are proposed and compared in 
I21i r3 1. Among these propagators, those based on the REGULAR con- 
straint do the most pruning and are often fastest. Finally, Bartak has 
proposed a similar intersection encoding for propagating a sliding 
scheduling constraint [ 1 1 We have shown that this method is more 
general and can be used for arbitrary SLIDE constraints. 

9 CONCLUSIONS 

We have studied the CardPath constraint. This slides a constraint 
down a sequence of variables. We considered SLIDE a special case of 
CardPath in which the slid constraint holds at every position. We 
demonstrated that this special case can encode many global sequenc- 
ing constraints including AmongSeq, CardPath, REGULAR in a 
simple way. SLIDE can therefore serve as a "general-purpose" con- 
straint for decomposing a wide range of global constraints, facilitat- 
ing their integration into constraint toolkits. We proved that enforc- 
ing GAC on Slide is NP-hard in general. Nevertheless, we identi- 



fied several useful and common cases where it is polynomial. For 
instance, when the constraint being slid overlaps on just one vari- 
able or is monotone, decomposition does not hinder propagation. 
Dynamic programming or a variation of the dual encoding can be 
used to propagate SLIDE when the constraint being slid overlaps on 
more than one variable and is not monotone. Unlike the previous 
proposed propagator for CardPath, this achieves GAC. Our exper- 
iments demonstrated that using SLIDE to encode constraints can be 
as efficient and effective as specialised propagators. There are many 
directions for future work. One promising direction is to use binary 
decision diagrams to store the supports for the constraints being slid 
when they have many satisfying tuples. We believe this could im- 
prove the efficiency of our propagator in many cases. 
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